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Abstract 

From data collected around the resonance by the OPAL detector at LEP, a 
sample of Bg decays was obtained using D'-^"^ combinations, where the D~ was fully 
reconstructed in the ij^vr", K*'^K~ and KgK~ decay channels or partially reconstructed 
in the (/>£~i>(X) decay channel. These events were used to study Bg oscillation. The 
flavor (b or b) at decay was determined from the lepton charge while the flavor at 
production was determined from a combination of techniques. The expected sensitivity 
of the experiment is 4.1 ps~^. The experiment was not able to resolve the oscillatory 
behavior, and we deduced that the Bg oscillation frequency Amg > 1.0 ps~^ at the 95% 
confidence level. 
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1 Introduction 

The phenomenon of B — B mixing is well established. In the case of the Bj system, the mass 
difference, Amd, between the two mass eigenstates has been measured rather precisely |jl|. 
This mass difference gives the oscillation frequency between Bd and Bd. Although these 
measurements can be used to gain information on the CKM matrix element Vtd-, this is 
hampered by large theoretical uncertainties on both the meson decay constant, /b^, and the 
QCD bag model vacuum insertion parameter, B-^^ 0. This difficulty may be overcome if 
the Bg oscillation frequency, Amg, is also measured. In this case, the CKM information can 
be extracted via the relation: 

Ams _ me, \Vts? fl^B^, 
Amd ruB^ |Vtdp Ib^Bb^ ' 

where ttlb^ and rriB^ are the Bg and Bd masses, as the ratio of decay constants for Bd and 
Bg mesons is much better known than the absolute values 0, |^. Information on |Vi;d| could 
then be extracted by inserting \Vts\, which is relatively well known [|l|. 

Arris is predicted to be many times larger than Amd [0; 0] and current lower limits support 
this theoretical predictions. A large Amg value leads to rapid oscillation thus presenting 
experimental difficulties, which have prevented its measurement to date. The most restrictive 
of the published limits [H, H, |^ indicates that Amg > 9.6 ps^^ at the 95% confidence level 0, 
while the best limit from OPAL gives Amg > 5.2 ps^^ at the 95% confidence level 

This paper describes an investigation of Amg using a sample enriched in Bg by recon- 
structing combinations^. In OPAL this technique is expected to achieve a sensitivity 
similar to that achieved by the inclusive technique 0], since the better decay time resolution 
and higher purity are offset by the lower statistics of an exclusive analysis. 

^Throughout this paper charge conjugate modes are imphed. 
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2 Analysis overview 



The oscillation frequency of Bg mesons was studied using exclusive decays of Bg mesons into 
D~£'^ combinations. Bg mesons were reconstructed in the following four D~ decay channels 
as described in 0. 

Bs ^ Dg- 1 V 

U 0rz/(X), ^ K+K- 



The selection procedure of the event sample followed closely that of and is briefly described 
in Section HI with an emphasis on the changes made to suit the purpose of an oscillation 



measurement. The background to the Bg signal is described in Section 

For each candidate we assigned a probability that it has mixed, i.e., its flavors (b or b) at 
production and at decay differ. This probability was derived from the decay and production 
flavor tags, and we refer to it as a mixing tag (Section ^. 

In order to assign a likelihood of a candidate at a given Amg value we need, in addition 
to the mixing tag, to reconstruct its decay time. Since the oscillation measurement is highly 
sensitive to the decay time, we did not assume a fixed Gaussian resolution on the decay 
time. We determined an event-by-event probability distribution for the decay time, which 
was derived from a Gaussian probability distribution for the decay length (Section and 
from a non-Gaussian probability distribution for the B candidate momentum (Section p.2| ). 

In order to extract a lower limit on the oscillation frequency, Amg, and to facilitate 
combination with other analyses we used the amplitude fit method |p. This method fits, 
for each value of Amg checked, a continuous parameter A which measures the size of the 
component in the data oscillating at that particular value of Amg. At the true Amg, the fitted 
value of A should be consistent with one, while far below the true Amg, the expectation 
value for A is zero (see [^j for additional details). Therefore values of Amg where A is 
below one and inconsistent with one will be excluded. The likelihood function and the fit 
results are described in Sections |^ and |^. The systematic effects of the uncertainties on 
all the parameters used in the amplitude fit were estimated by repeating the amplitude fit 
with those parameters varied by one sigma (Section Several checks of the method are 
described in Section [TO. Finally our results, and the results of combining this measurement 



with the previous OPAL measurement are summarized in Section 11 



3 Hadronic event selection and simulation 



We used data collected by the OPAL detector [|10[ at LEP between 1991 and 1995 running 
at center-of-mass energies in the vicinity of the Z° peak with an operational silicon detector. 
Hadronic Z° decays were selected using the number of tracks and the visible energy in each 
event as in [|1T|. This selection yielded 4.3 million hadronic events. In each event, tracks 
and electromagnetic clusters not associated to a track were combined into jets, using the 
JADE algorithm with the EO recombination scheme. Within this algorithm jets are defined 
by|/,,t = 0.04 ig. 

Monte Carlo samples of inclusive hadronic Z° decays and of the specific decay modes 
of interest were used to check the selection procedure, mix tagging and fitting procedure. 
These simulated event samples included: 
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Samples of the four signal decay channels. 



Hadronic T? decay samples, used to check the selection efficiencies and mix tagging of 
Z° — > qq decays, where q is a light quark (u, d, s or c). 

Z*^ — >■ bb decay samples, used to check the selection efficiencies and the mixing tag of 
other background decays, such as partially reconstructed signal decays (Section |6.3.1| ). 

Z° decay samples containing the specific decays Bd D+D^, Bu — > D+D*^, Bd — > 
D+D-**, Bu ^ D+DO", Bu ^ D+DOtt, Bd ^ D-K^^+i/ and B^ ^ D-K+£+z/, used to 
check the selection efficiencies of D~£"'" background decays (Section [4.2. 



These samples were produced using the JETSET 7.4 parton shower Monte Carlo gener- 
ator [jl3| with the fragmentation function of Peterson et al. for heavy quarks, and then 



passed through the full OPAL detector simulation package [15 



4 Candidate selection 

Three tracks were combined to form a D~ candidate and a lepton (either e or /i) was added 
to form a Bg candidate. The four tracks were required to be in the same jet. 

The event selection and decay length reconstruction for this analysis follow closely those 
of 0. Since the decay time resolution is crucial in this analysis, an additional requirement 
was made, demanding that the prompt lepton track (that is the lepton directly from the 
B decay) had at least one associated hit in the silicon microvertex detector. The photon 



conversion rejection has been updated to use a neural network ||T6|. The event selection and 
reconstruction are outlined briefly below: 

Standard track quality cuts il^ were applied. Electrons were identifled using a neural 



network and a photon conversion rejection cut; muons were identifled by associating 



central detector tracks with track segments in the muon detectors and requiring a position 



match in two orthogonal coordinates [|l^. For the other reconstructed particles the proba- 
bility that the observed rate of energy loss due to ionisation (d-E/dx) is consistent with the 
assumed particle hypothesis was required to be greater than 1%. 

Additional channel dependent cuts included: momentum cuts, further di?/dx cuts, invari- 
ant mass cuts on reconstructed intermediate particles, including a loose cut on the invariant 
mass of the visible Bg decay products, helicity angle cuts and a cut on the angle between 
the D~ candidate and the prompt lepton candidate. See [0 for details. 

In the Kg K~ channel the mass of the two tracks forming the Kg candidate was constrained 
to the known Kg mass Further constraints were applied to the D~ and Kg, in which the 
directions of the vectors between their production and decay points were constrained to the 
reconstructed momentum vectors. The lepton minimum momentum cut in this channel was 
5GeV. 

Three vertices were reconstructed in the x-y planeQ: the e^e" interaction vertex, the Bg 
decay vertex and the D~ decay vertex. The e^e~ interaction vertex was measured using 
tracks with a technique that follows any signiflcant shifts in the e^e' interaction vertex 
position during a LEP flll [BO]. The D^T decay vertex was fltted in the r-0 plane using all 



^The right-handed coordinate is defined such that the z-axis follows the electron beam direction and the 
x-y plane is perpendicular to it with the x-axis lying approximately horizontally. The polar angle 9 is defined 
relative to the -|-z-axis, and the azimuthal angle is defined relative to the -l-x-axis. 
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the candidate tracks. The Bg decay vertex was formed by extrapolating the candidate D~ 
momentum vector from its decay vertex to the intersection with the lepton track. 

The D~ decay length is the distance between these two decay vertices. The Bg decay 
length was found by a fit between the e^e~ interaction vertex and the reconstructed Bg 
decay vertex using the direction of the candidate D~£"'" momentum vector as a constraint. 
The two-dimensional projection of the Bg decay length was converted into three dimensions 
using the polar angle that was reconstructed from the momentum of the D~£"'". Typical 
reconstructed decay length errors range from about 0.35mm for the K"'"K~7r~, and 0£~z/ 
channels, to about twice this level for the K^K" channel. In channels where the D~ was 
fully reconstructed the of the D~ decay vertex fit was required to be less than 10 (for one 
degree of freedom). Finally, the reconstructed decay length error of the Bg candidate was 
required to be less than 0.2 cm. 



4.1 Results of Dg selection 

The invariant mass distribution obtained in each of the D~ decay channels is shown in Fig- 
ure m. Each invariant mass distribution was fitted to a Gaussian distribution describing the 
signal and a linear parameterization for the combinatorial background. The mean of the 
Gaussian distribution was fixed to the nominal D~ mass, 1968.5 MeV [||], for the hadronic 
channels and to the nominal mass, 1019.413 MeV [|ll, for the semileptonic channel. In 
the K+K"??" distributions, a second Gaussian distribution was used to parameterize con- 
tributions from the Cabibbo suppressed decay D~ — K+K"??^. The mean of the second 
Gaussian distribution was fixed to the nominal D~ mass, 1869.3 MeV and the width was 
constrained to be the same as that of the D~ peak. The combinatorial background in the 



semileptonic channel was refitted to account for the kinematical threshold as in |]2T|. The 
choice of the background parameterization was found to have a negligible effect on the fitted 
amplitude. For each channel, the fitted width was consistent with the expected detector res- 
olution. The contamination due to the D~ — >■ K"'"K~7r~ decays was estimated from simulated 
events as explained in Section |6.3.1| . The results of these fits are summarized in Table |I|. 



No significant peaks were observed in the mass distributions for same-sign D~£~ combi- 
nations in the fully reconstructed decay channels K*'^K^, 07r^ and KgK^. 

The signal and sideband regions are defined in Table ^ and shown in Figure |1]. The D~£"'" 
combinations selected for the oscillation fit were from the signal regions. 244 such D~£"'" 
candidates were observed. The candidates selected in the mass sideband regions were used 
to estimate the lifetime characteristics of the combinatorial background. 199 such sideband 
candidates were observed. The regions were chosen to represent the decay time distribution 
of background under the signal as in |^. 



4.2 Background to the Bg — > signal 

Potential sources of non-combinatorial background to the Bg signal considered here include 
decays of other B hadrons that can yield a D~£"^ final state or other final states that were 
misidentified as a D~ meson. Other sources are a D~ combined with a hadron that has been 
misidentified as a lepton, and random associations of a D~ with a genuine lepton. Finally, 
there is combinatorial background from misreconstructed D~ mesons in the hadronic chan- 
nels and misreconstructed mesons in the semileptonic channel. The various background 
sources, and the calculation of their contributions relative to that of the signal, are discussed 
below. 
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OPAL 




1.6 2 2.4 0.98 1.03 1.08 

m(K*K) (GeV) m(KK) (GeV) 

Figure 1: Invariant mass distributions from the four D~ decay channels. In each plot, the 
result of the fit to the signal and possible satellite peaks, as described in the text, is overlaid 
as a solid line. The bold arrow to the left shows the signal region and the lighter arrow to 
the right shows the sideband region. 
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Decay 
channel 


Candi- 
dates 


Comb, 
fraction 


Signal 
region (MeV) 


Sideband 
region (MeV) 


Estimated 
signal 


K*K 

(f>e 


125 
54 
24 
41 


0.455±0.026 
0.277±0.027 
0.467±0.092 
0.243±0.039 


1918.5 - 2022.1 
1929.5 - 2007.5 
1918.5 - 2091.3 
1011.4 - 1027.4 


2022.1 - 2168.5 
2022.1 - 2168.5 

2091.3 - 2168.5 

1027.4 - 1079.4 


53.8±5.0 
30.9±2.7 
10.1±2.3 
21.4±2.8 


Total 


244 


0.386±0.018 






116±10 



Table 1: Results of the mass fits to all the signal channels. The second column shows the 
number of events selected in the signal region as defined in the text. The fitted fraction of 
combinatorial background for events selected in the signal region is given in the third column, 
together with the associated statistical error. The fourth column gives the mass range for the 
signal region, which corresponds to about twice the fitted width around the nominal D~ mass 
for the hadronic channels, and to about twice the width around the nominal (f) mass for the 
semileptonic channel. The fifth column gives the mass range for the sideband region. The 
last column gives the estimated number of correctly reconstructed Bg ^s^~^ signal events 
in the signal region, where both combinatorial and other background discussed in [j.^ were 
subtracted. The errors are from the uncertainties in the combinatorial and non- combinatorial 
background subtraction. 



4.2.1 T>-1+ Background 

The event sample includes properly reconstructed D~£"^ combinations that do not arise from 
Bs decay. Two decay modes of Bj and Bu mesons were considered: 

(a) B(u/d) DjrD(X), D £^z/(X) (where D is any non-strange charm meson). 

(b) B(u/d) DjrK£+z/(X), where K also represents excited kaons. 

Note that in decay mode (a) a negative lepton would indicate a B^ or Bj meson, whereas in 
signal decay and in decay mode (b) a negative lepton would indicated a Bg, Bu or Bj. 

Mode (a) includes two body decays B(u/d) D^rD for which a branching ratio measure- 
ment exists, and three body decays B(u/d) DfT-DX for which no branching ratio measure- 
ment exists. The measurement Br(B(u/d) D^rX) = 0.100 ± 0.025 IQ provides an upper 
limit on the sum of these modes. To estimate this mode's contribution to the D^f"*" sam- 
ple, the semileptonic branching ratios for the different non-strange D hadrons are weighted 
according to their abundance in B(u/d) — > D~D decays. The efficiency for this analysis to 
select D~£"'" combinations from these modes was taken from Monte Carlo and included in 
the calculation. The possibile contribution from analogous decay modes of b-baryons into 
was included in the calculation although none have been observed to date. These 
modes account for 0.143 ± 0.050 of the selected D~£+ combinations, and 0.619 ± 0.045 of 
this background comes from B^ decays. 

Mode (b), B(u/d) D~K£+z/(X), has not been observed and only an upper limit of 0.009 
(90%CL) m exists. However, a branching ratio for this mode can be calculated from the 
fraction of 'lower vertex' (D^T not produced from the virtual W) B Dg^X ||l|| combined 
with the inclusive B semileptonic decay fraction. This yields Br(B D~K£"*"i/(X)) = 
0.0018 ±0.0009 which is consistent with the above upper limit as well as with the theoretical 
upper limit [E^. Monte Carlo events were used to determine the selection efficiencies for 
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these background modes relative to that of the signal mode. Analogous baryonic modes were 
included in this calculation as well. These modes account for 0.065 ± 0.035 of the selected 
D~^+ combinations, and 0.470 ± 0.038 of this background comes from decays. 



4.2.2 Other background 

The background from genuine D~ particles that were combined with a hadron that was 
misidentified as a lepton can be estimated from the invariant mass spectrum of combinations 
of same-sign charm candidate and lepton candidate pairs. Assuming that misidentified 
hadrons are equally likely in both charges, the same number of D~-|-fake lepton should exist 
with the correct charge correlation. For each channel in which the charm hadron is fully 
reconstructed, no significant excess of same sign D~ signal exists. This is in agreement with 
what has been found in a related analysis that has greater statistical significance |2^. This 
background source was therefore neglected. 

In the channel </)£+z/(X), where the charm hadron was partially reconstructed from 

a semileptonic decay channel, there was additional background to consider. This background 
includes the accidental combination of a 0, produced in fragmentation, with two leptons 
that arise from either B — > D£X, D — > £X or B — >■ J/-?/' decays and candidates from hadrons 
misidentified as leptons. In it was estimated that the fraction of candidates in the (p 
signal region that arise from this background particular to the 0£~''z/(X) channel is 

/other = 0.135 ± 0.057, an estimate used here as well. 

The non-combinatorial background sources mentioned above were expected to contribute 
a total of 35 ± 9 events to the D"^"*" signal. The background subtracted number of D~£"'" 
signal candidates was therefore N{Bs — > D~£+z/(X)) = 116 ± 10, as given in Table |1|. 



5 Proper decay time reconstruction 

The true Bg proper decay time, t, is derived using the relation: 

hrne ' 



(2) 



where /true, Ptrue and m-Q are the Bg candidate's true decay length, true momentum and 
nominal mass respectively. From the measured decay length we derive a Gaussian probability 
distribution for Ztruo, as described in Section |0| . We also derive a non-Gaussian probability 



distribution for ptrue, as described in Section 5.2 



5.1 Decay length estimation 

The Bg candidate's decay length is reconstructed as described in Section ^ Using simulated 
events it was found that the decay length reconstruction is biased and that the decay length 
errors reconstructed by this method were overly optimistic by a factor of about 1.4. On 
average the reconstructed decay length was bigger than the true decay length by 24 ±12 ixm. 
We corrected for this bias, which is 6% of the average decay length resolution, and less than 
1% of the average decay length. The distribution of the reconstructed decay length errors 
in the data and in simulated events is similar, as shown in Figure |^. Using simulated signal 
events we fitted the ratio between the correct decay length error, a\ and the reconstructed 
decay length error, 0"'^^^^, as a linear function of crl^con- Simulated signal events from all 
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Figure 2: Distribution of reconstructed decay length errors for selected events. The points 
with the error bars represent the data while the histogram represents the simulated events. 

four decay channels were used, and the dependance of a' on cr^g^^^ was similar in all signal 
channels. We used this function to correct the decay length error in the likelihood function 
calculation, and used the fitted uncertainty on this function as a systematic error. 

5.2 Momentum estimation 

Since the prompt neutrino produced in the candidate's decay, and in some cases additional 
decay products, are not reconstructed, there is no direct measurement of the candidate's true 
momentum ptrue- The binned probability distribution of the candidate's true momentum, 
B{ptrue), is estimated on an event-by-event basis using a probability distribution based on 
the reconstructed candidate {Bi, Section [5.2. 1|) and a probability distribution based on 



the recoil to the candidate, i.e. the other tracks and clusters in the event {B2, Section [5.2.2| ). 
The two probability distributions were then used to calculate B using: 

. _ ^l(Ptruc) • -g2(Ptrue) . . 

where n is the number of momentum bins. 

5.2.1 Candidate based momentum distribution {B\) 

We calculate a probability distribution for the Bg candidate's true energy, i?B, using the 
reconstructed invariant mass, mo^, and energy, -Ed^, of the D~ lepton combination as exper- 



10 



imental inputs, following the method presented in [23 . 

A Bayesian approach is used for which an a priori knowledge of the Bg candidate's en- 
ergy spectrum is required. This a priori spectrum, P{E^), was derived from Monte Carlo. 
Applying two body decay kinematics, the observable energy, E^i, is given in the laboratory 
frame by: 

^d£ = 77^(S + /?bAcos^*) , (4) 
IrriB 

where 6^ is the angle between the D£ flight direction and the boost vector in the Bg candi- 
date's rest frame, S = m| -|- m^^ , A = m| — m^^ , Pb and 7b are the boost parameters of 
the Bs candidate in the laboratory frame. 

The distribution in cos6'g is uniform (because the B meson is a pseudoscalar particle), 
therefore E^i is distributed uniformly between 2^(^ ~ Pb^) and 2^(^ + Pb^)- 

We then used the fact that E^ is independent of the Di invariant mass to get P{Eb, rriDi) = 
P{E-q) ■ PirriDi), together with Bayes theorem to obtain the formula: 

p(T^ \-p ^_ P{Em\EB,mm)P{EB,mm) _ P{EY^i\EB,moe)P{EB) 

[ Bl ^^^'^^^)- jp^E^^\E'^,m^e)PiE'^,m^e)dE'^~ JPiE^e\EB,mne^^^ ^ ^ 

The momentum probability density, Bi, is then derived from the energy probability 
density. Using simulated signal decays, it was found that on average the expectation value 
of Bl was smaller than the true momentum by 0.24 ± 0.06 GeV. We corrected for this bias, 
which is less than 1% of the average candidate momentum. 

5.2.2 Recoil based momentum distribution (B2) 

Another way of obtaining a good estimate of the Bg candidate's momentum is to use our 
knowledge of the total center of mass energy Ecm, which is twice the LEP beam energy. We 
calculate the Bg candidate's energy using this constraint and the recoil mass of the rest of 
the event, using the relation: 

^cm^ + m| - mg,, 
Eb = , (6) 

where mrec is the recoil mass calculated using all tracks and unassociated electromagnetic 
clusters in the event except the reconstructed Bg decay products, and tub is the nominal Bg 
mass. In this calculation all tracks were assigned the pion mass and all neutral clusters were 



taken as massless. The candidate's momentum is then given by y E"^ — m|. 

Monte Carlo studies have shown that the accuracy of this estimate can be improved by 
rescaling mrec according to the visible energy, calculated using all tracks and unassociated 
electromagnetic clusters in the event, and Ecm (see Figure |^). 

Figure ^ shows that the difference between the corrected reconstructed momentum and 
the true simulated momentum (ptme — Precon) is well described by a Gaussian distribution 
whose center is at 0.31 ± 0.05 GeV. We corrected for this bias, which is less than 1% of the 
average candidate momentum. Therefore B2 was chosen as a Gaussian distribution around 
the reconstructed momentum minus the bias (precon — 0.31 GeV) with a width of 2.88 GeV, 
i.e. the fitted Gaussian width of ptrue — Precon in simulated signal events. 

5.2.3 Results of momentum estimation 

The width (RMS) of the Bi distribution varies greatly between events, the average width on 
simulated signal events is 3.7 GeV. As stated above the B2 distribution has a single width 
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difference between constrained and true P^^^^^ (GeV) 

Figure 3: Comparison of errors in the momentum estimate B2 before (upper) and after 
(lower) the visible energy correction of m-cec on simulated signal events. The points with the 
error bars represent the simulated events while the solid line shows a Gaussian fit to the 
points. 

of 2.88 GeV for all events. The average width of the combined distribution B on simulated 
signal events is 2.30 GeV. 

The small individual biases on Bi and B2 were corrected before combining them according 
to Equation |^. It was verified on the simulated signal events that after correcting for both 
biases, B is indeed a reasonable representation of the true probability distribution for the 
candidate's momentum. 

5.3 Results of proper decay time estimation 

The distribution of the true proper decay time is estimated by combining the decay length 
estimate (|5.1|) and the momentum estimation (|5.2| ) according to Equation 0, after correcting 
for their small biases. Figure ^ shows that for 68% ±5% of simulated signal events the differ- 
ence between the expectation value of the reconstructed true proper decay time distribution, 
texp, and the true decay time, t, is well characterized by a Gaussian distribution of width 
0.175 ± 0.011 ps. This fit is shown for information only, and was not used in the oscillation 
fit likelihood. 

We classified simulated signal events according to the RMS of their reconstructed true 
proper decay time distribution, a*. We found that in events with low a* the expectation 
value of the reconstructed true proper decay time distribution, tgxp, tends to be smaller than 
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Figure 4: Difference between the expectation value of the estimated probability distribution 
of the true proper decay time, texp; and the true proper decay time, t. The points with the 
error bars represent the simulated events while the solid line shows a double Gaussian fit to 
the points. 68% of the events lie within the narrow Gaussian distribution whose width is 
0.175 ps. The width of the wide Gaussian distribution is 0.52 ps. 

the true proper decay time, t. In events with high a*, tcxp tends to be bigger than t. This 
residual bias is about 2% of the average proper decay time. We fitted the proper decay 
time reconstruction bias as a hnear function of cr* (^bias'^'^* + ^bias)' ^^"^ treated the fitted 
uncertainties as sources of systematic uncertainty. 



6 Mixing tag 

The mixed and unmixed Bg decays were distinguished by determining the b flavor of the Bg 
(whether it contains a b or b quark) both at production and at decay. The b flavor at decay 
was inferred from the charge of the prompt lepton in the Dg"^+ combination. The initial 
b flavor was tagged by a combination of the charge of a lepton in the hemisphere opposite 
the Bg candidate, the charge of a fragmentation kaon in the candidate hemisphere, and jet 
charge measures from both the candidate hemisphere and the opposite hemisphere. The 
available tags in each hemisphere were combined into a measure of the probability that the 
candidate was produced as a Bg and then the probabilities from the two hemispheres were 
combined into a single probability. It was verified that any unwanted correlations between 
the flavor tags of the two hemispheres were negligible. The mixing probability is derived 
from the production flavor probability and the decay flavor. 
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6.1 The Bg candidate hemisphere 

The Bg production flavor was measured in the Bg candidate's hemisphere by the jet charge 
and, where available, the charge of a kaon from the fragmentation process. 
The jet charge of the jet containing the Bg candidate was calculated as 

1 " 

Qsame = " E " (p'T , (7) 

where p\ is the longitudinal component of the momentum of particle i with respect to the 
jet axis, Qi is the electric charge (±1) of particle i. The sum is over all tracks in the jet 
excluding the Bg decay products, since the latter contain no information on whether the 
candidate meson was produced as a Bg or Bg and would only dilute the information from 



the fragmentation tracks. The optimal value of k was found to be 0.4 as in [gj . 

The fragmentation kaon tag is an attempt to identify the kaon containing the s quark 
that was produced in the fragmentation process in association with the s quark which is part 
of the Bg. This kaon was selected as follows: 

• The probability that the measured dE/dx of the track is consistent with the kaon 
hypothesis is greater than 1%. 

• The measured dE/dx of the track is lower than the expected value for pions of that 
momentum by at least one standard deviation of the measurement. 

• The measured dE/dx of the track is higher than the expected value for protons of that 
momentum by at least one standard deviation. 

• The track is not identified as a lepton (as in Section ^). 

• The distance of closest approach in r0 of the track to the e^e^ interaction vertex is 
smaller than 2 mm. 

When two tracks with the same reconstructed charge satisfied these requirements, the event 
was tagged using that charge. When the two tracks' charges were different, or when three 
or more tracks satisfied those requirements, the event was not given a fragmentation kaon 
tag. The latter scenario is limited to less than 5% of the tagged events. 

When no fragmentation kaon was tagged, the jet charge was converted to a probability 
using the Bayesian formula: 

P (Bg I Qsamc) 



P (Qsamc 


Bg) 


P(Qsamc 


Bg) ~l~ P(Qsame 


Bg) 



where P(Bg) is the probability of the candidate being a Bg and not a Bg, and P{Qsa.me\fiO'Vor) 
is a Gaussian probability density describing the Qgame distribution conditioned by the candi- 
date's true production flavor, as obtained from a fit to signal Monte Carlo. This formula uses 
the fact that the a priori probabilities of both flavors are one half. The separation between 
the two competing hypotheses is shown in Figure 

When an additional fragmentation tag was found, the jet charge and the additional tag 
were converted to a probability using the following Bayesian formulae: 
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P (BslQsame, T = Bg 



P(Qsamc; Bs|Bs) -|- P(Qsame; Bs|Bs) 



P (Qsamc 1 


tag) ■ P(tag) 


P('5same 


tag, B,) ■ P(tag) + P(Qsame 
P(Qsame)-^ Bs|Bs) 


mistag, Bs) ■ P(mistag) 



V> fn \n T — n\ — Wsame, J- — -Dsl-Dsj /^n 

-•^ I, same 5 -DslJDgJ ~r -T ^^li^/game; -Dsl-Ds/ 



P (Qsame 


mistag) ■ P (mistag) 


P(Qsame 


mistag, Bs) ■ P (mistag) + P((5same 


tag, Bs)-P(tag) 



where T is the flavor indicated by the kaon tag, and P((5same|mltag5 ^s) are two Gaussian 
probabihty distributions fitted on simulated Bg decays for the case when T indicates the 
correct flavor and for the case when the wrong flavor is indicated, and P(Qsamelmfstag5 ^s) = 
P(— Qsamelmltag) Bs). Again usc was made of the fact that the apriori probabilities of both 
flavors are one half. A comparison of the two competing hypotheses is shown in Figure §3. 

The Qsame distribution is not necessarily charge symmetric because of detector effects 
causing differences in the reconstruction of positively and negatively charged tracks. These 
effects are caused by the material in the detector and the Lorentz angle in the jet chamber. 
They were removed by subtracting an offset from the Qsame value before using it to tag 
the candidate's production flavor and before parameterizing P((5same|miitag) Bs). The small 
Qsame offsct was determined from simulated signal events, since no pure sample of fully 
reconstructed signal decays is available from the data. This procedure gains support from the 
agreement between the Qopp offset values calculated from simulation and data in Section |0[ 
After subtracting the offset, the simulated Qsame distribution is charge symmetric. The Qsame 
offset was found to be 0.006 ± 0.004, where the error is from limited Monte Carlo statistics. 

6.2 The jet opposite the Bg candidate 

Flavor ant icor relation between the two hemispheres allows the use of the b flavor in the 
hemisphere opposite the candidate to tag the candidate's production flavor. The b flavor in 
that hemisphere was tagged using the jet charge of the highest energy jet it contains, and, 
where available, the charge of a track identifled as a lepton from semileptonic b decay. 

The jet charge in the highest energy jet opposite the Bg candidate, Qopp, was calculated 
in the similar way to Qsame, except that here the sum included all the particles in the jet and 



the optimal value of k was found to be 0.5 as in [^. This value of k optimizes the weight 
given to the fragmentation tracks' charges relative to the weight given to the decay tracks' 
charges. 

A lepton in the opposite hemisphere was selected as follows: 

• The track is identified lepton as in Section H. 

• Momentum greater than 2 GeV. 

• Transverse momentum greater than 0.8 GeV with respect to the jet axis. 

• It must not be identified as arising from photon conversion. 

When two tracks with the same reconstructed charge satisfied these requirements, the event 
was tagged using that charge. When the two tracks' charges were different, or when three 
or more tracks satisfied those requirements, the event was not given an opposite lepton tag. 
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Figure 5: Distribution o/Qsame forBs candidates in simulated events. Points with error bars 
show the number of events, the solid lines show the Gaussian fits used in the analysis and 
the broken lines show the relevant competing hypothesis. The top plot shows events without 
a fragmentation kaon tag, where the competing hypothesis is the opposite candidate flavor. 
The lower plot shows events where the fragmentation kaon tag indicates the correct candidate 
flavor, where the competing hypothesis is the opposite candidate flavor and fragmentation 
kaon mistag. Note that the competing hypothesis has a lower overall probability since the 
fragmentation kaon tag's purity is more than 50%. 



The jet charge and the lepton tag of the opposite hemisphere were converted to a probabil- 
ity using the same method as in Section |6.1| for the hemisphere containing the Bg candidate. 
A comparison of the two competing hypotheses is shown in Figure ^. 

As described in Section ^]T] for Qsame, the Qopp distribution is not charge symmetric. 
The Qopp offset was determined using a large sample of b tagged inclusive lepton events 
selected from data. The resulting value of the Qopp offset agrees well with the values derived 
from simulated signal events, from a large simulated sample of b tagged inclusive lepton 
events, and from [Q. After subtracting the offset, the simulated Qopp distribution is charge 



symmetric. The Qopp offset was found to be 0.0138 ± 0.0020, where the error is from the 
limited statistics of the selected data sample. 
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Figure 6: Distribution o/Qopp for Bg candidates in simulated events. Points with error bars 
show the number of events, the solid lines show the Gaussian fits used in the analysis and the 
broken lines show the relevant competing hypothesis. The top plot shows events without an 
opposite lepton tag, where the competing hypothesis is the other candidate flavor. The lower 
plot shows events where the opposite lepton tag indicates the correct candidate flavor, where 
the competing hypothesis is the other candidate flavor and opposite lepton mistag. Note that 
the competing hypothesis has a lower overall probability since the opposite lepton tag's purity 
is more than 50%. 



6.3 Mixing tag results 

The procedure described above attempts to assess the probabihty that an event underwent 
mixing. Being based on Gaussian approximations of the jet charge distributions, this raw 
mixing tag, x, is therefore only an approximation of the true mixing probabihty. A calcula- 
tion of the events likelihood demands that we calibrate the mixing tag, M , by quantifying 
its deviation from a true probability. 

The average mixing tag for simulated signal events, in which half the decays were mixed, 
was found to be consistent with one half. Furthermore the difference between this average 
and the true average (0.5) is much smaller than the typical systematic uncertainties on the 
mixing tag, and is neglected in this analysis. The distribution of the mixing tag in simulated 
signal events is shown in Figure 0. 
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Figure 7: Distribution of the raw mixing tag in simulated signal events. The light histogram 
shows the distribution of the raw mixing tag, the dark histogram shows the contribution of 
mixed events. 

The deviation of the mixing tag from the a posteriori probabihty is parameterized with 
the parameters a^ix and P^nix (used in the oscillation fit Section which quantify the 
deviations at x = 0.75 and x = 1.0 respectively. The resulting fit is shown in Figure ^ and 
the fitted uncertainty is treated as a systematic error. 

6.3.1 Mixing tag behavior in the combinatorial background 

To calculate the likelihood of an event originating from combinatorial background, we need to 
know the behavior of the mixing tag in combinatorial background events. Two effects which 
influence this behavior were found using Monte Carlo: a different mixing tag distribution, 
and oscillations of the background. These effects were identified in particular subsamples of 
the combinatorial background. The definition, abundance and behavior of those subsamples 
follows. 

Combinatorial background events which originated in bb events can exhibit oscillatory 
behavior if the decaying meson is a Bd or a Bg. Two rates of oscillation were found in 
simulated combinatorial background events: 

Oscillation with the same rate as the simulated signal oscillation rate arises primarily 
when the decay is truly through a Bs and only one of the D~ decay products is misidentified. 
The proper decay time reconstruction for these simulated events was less accurate than for 
signal events, the width (RMS) of the difference between the true time and the expectation 
value of the reconstructed time distribution is bigger by 25 ±7% for those events. The mixing 
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Figure 8: Deviation of mixing tag from a true probability. The points with the error bars 
show simulated events, the solid line is a fit to these points as described in the text and the 
dotted line is M = x, representing a true probability requiring no correction. The error bars 
include systematic errors described in Section Both charge conjugated cases were added 
to maximize the statistical significance. 

tag for these simulated events was essentially the same as for simulated signal events. 

It was found that oscillation with the same rate as the Bd oscillation rate arises primarily 
when the decay products of a D+ or D° meson produced from a Bd were reconstructed 
as a D~ meson, and at most one of the D meson decay products was misidentified. A 
typical decay chain for this channel is Bd D^*H~^i'(K), which is similar to the signal 
decay Bg — > Ds£"'"z/(X). The mixing tag performance and decay time reconstruction for these 
simulated events were consistent with that for simulated signal events. 

Monte Carlo studies have shown that the fraction, C, of combinatorial background which 
oscillates depends on whether the channel contains a 0, both for Bg and Bd fractions. There- 
fore we use four parameters in the fit, Cq, Cq° Cf and Cf° ^, that give the Bg and Bd 
fractions for (j) and No-(/) channels. 

The distribution of the mixing tag for the non-oscillating combinatorial background is 
different from that for signal, hence the mixing tag not only indicates mixing but also 
contributes information as to whether the event is a signal event. Our use of this information 
is described in Section |6.3.2| . 

Biases were found mainly on non-bb events. While uds events tend to be tagged as mixed, 
heavier flavor events have a statistically significant tendency to be tagged as unmixed. These 
tendencies partially cancel out, leaving an overall bias of the mixing tag on non-oscillating 
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combinatorial background events that was found to be B^°^ = —0.030 ± 0.009. 
6.3.2 Mixing tag behavior in the D~£+ background 

The mixing tag for simulated D~^+ background events of type (a) (as defined in Section ^. 2. 1| ) 
was essentially opposite that for simulated signal events, as expected from the decay chain. 
In addition, the distribution of the mixing tag for D~£+ background from Bu decay was 
found to be different from that for signal, while for D~£'^ background from Bd decay the 
distribution of the mixing tag was consistent with that for signal. Hence the mixing tag 
not only indicates mixing, but also contributes information as to whether the event is from 
the signal, from the combinatorial background (Section |6.3.1| ), or from the Bu 
background. We derived a probability, /™'^, given only the mixing tag, that an event is 
from the combinatorial background, and similarly a probability, ff^^^^ that it is from the 
Bu D^^^ background. Using simulated events we fitted f™^^ and ff^^^ as functions of the 
mixing tag, and the fitted uncertainties were used as a systematic uncertainties. 



7 Oscillation fit 

The likelihood, £, for observing a particular decay length, /j, of candidate i, and a particular 
mixing tag, Mj, may be parameterized in terms of the candidate's decay length error, cr-, its 
calculated Bg momentum spectrum (Section ^3), Bi, and a probability, /™^^^, that it arises 
from combinatorial background. /™^^^ is determined as a function of the observed invariant 
mass of this candidate from the fit to the invariant mass spectrum shown in Figure |1|. 

An event's likelihood is found by summing over all the possible event types (i.e. signal, 
the two D~£"'" background modes, the two oscillating combinatorial background modes, and 
regular combintorial background). For each event type we assign a probability that it is an 
event of this type, P(type), and the likelihood if the event is of that type, C^i'^'^: 

fcomb / /mass i'mixN i'physX /'type/? Tir i / 



event type 

(10) 

where f^'^'^^ is the estimated probability that the event arises from combinatorial background 
based on the invariant mass (through f^^^^^ and on the mixing tag (through /j™^), and ff^^^ 
is the estimated probability, based on the mixing tag, that the event is from the Bu — >■ D^^"*" 
background. 

The form of the likelihood function for signal events is given by the convolution of three 
terms: a term describing the probability of the mixing tag and the true decay length given the 
true momentum, the calculated momentum distribution, and a Gaussian resolution function 



with width equal to the decay length error (corrected as described in Section This can 
be expressed as: 

Cr''\luM,\alB;) = / d/true,i (11) 

/ dptfuc QiJ-i I ^true,i) ^i) ^j(Ptrue) "^signal (^true,j ; -^^i | Ptrue) ; 
J 

where the function ^ is a Gaussian function that describes the probability to observe a decay 
length, /j, given a true decay length Ztrue,i and the estimated measurement uncertainty a\. 
Bi{ptiue) is the probability of a particular Bg momentum. V is the probability for a given Bg 
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to decay at a distance /true,i from the e'^e interaction vertex with a mixing tag Mj. This 
function is given by: 

-P (J Mir, \ - exp(-t/rBj 

Ptrue r-B, 



1 + A- cos(tAms) 1 — ^ ■ COS (t Arris) ,^ 
^ ■ (1 - Mi) + ^ ■ M, 



where tb^ is the Bg hfetime, ^ is the fitted amphtude of the oscillation and the Bg 
candidate's true proper decay time, t, is given by Equation |^. 

For the D~£+ background events from Bj decay, the likelihood is similar except that the 
Bd lifetime was used to obtain: 

V (I ■ M- I « 1 - ""^ exp(-t/rBj 

' osc V''true,i 5 -"^-'j | f true J 



Ptrue Tb, 



d 



+ cos(tAm,) . _ ^ 1 - cos(tAm,) _ ^^^^ 

where for decay mode (a) we need to replace Mj with 1 — Mj. For the D~£"'" background 
events from Bu decay, the likelihood is simpler, and contains an exponential decay term with 
the Bu lifetime weighted by the candidate's probability to be unmixed. 

The combinatorial background was divided into several types according to oscillatory 
behavior, with most of the combinatorial background events being of the non-oscillating 
type. The function used to parameterize the reconstructed decay length distribution of this 
background is the sum of a positive and a negative exponential, convoluted with the same 
boost function as the signal and a Gaussian resolution function. This can be expressed as: 



/+ exp 
P (J, I + - /■+ \ _ ) '^g^bgPt™ 



;i - /bg) exp 



iili>0 

(14) 

if /, < . 



The fraction of background with positive lifetime, /|^, as well as the characteristic positive 
and negative lifetimes of the background, t^^ and r^^, were obtained from a fit to the sideband 
region. The resulting value and their uncertainties were used to constrain the background 
lifetime parameters in the oscillation fit. The background parameters were fitted separately 
for the hadronic and semileptonic D~ decay channels, as in |0. The lifetime behavior of 
the hadronic D~ decay channels' sidebands was best described by fitting only a positive 
exponential decay. The lifetime behavior of the semileptonic D~ decay channel's sideband 
was best described by fitting both exponential terms. 

For the semileptonic channels, the background which include a real not from a D~ is 
treated as combinatorial background. For the oscillating types of combinatorial background 



we used the following: for background oscillating at the Bd frequency we used Equation ^ 
while for background oscillating at the Bg frequency we used Equation |12[ 

8 Results of oscillation fit 

The results of the amplitude fit to the selected events are shown in Figure ^, including the 
systematic uncertainties (Section ||). An amplitude peak is evident at Amg = 6.0 ps~^, above 
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the experimental sensitivity, but nevertheless it seems inconsistent with an amplitude of zero 
with a significance of 2.35 sigma (including systematic uncertainties). The current combined 
world lower limit is Arris > 14.4 ps~^; this leads us to interpret this peak as a statistical 
fiuctuation. At low frequencies the fitted amplitude quickly rises above the A = line, and 
therefore after taking into account all systematic uncertainties described in Section ^, this 
analysis can only set a weak lower limit of Arris > 1.0 ps^^ at the 95% confidence level. 



9 Systematic uncertainties 

The systematic uncertainties on the Bg oscillation amplitude, (jsystematic^ ^^^^ calculated, using 
the prescription of p|, as: 

new nominal 

^systematic ^new ^nominal _|_ |^]^ ^nominal '\^_ fl5) 

^ ' ^nominal ' ^ ' 

where the superscript "nominal" refers to the amplitude value. A, and statistical uncertainty, 
a, obtained using the nominal values of the various parameters, and the "new" refers to the 
new values obtained when a single parameter is increased or decreased by its uncertainty and 
the fit is repeated. The systematics shown are an average of the effects of the increment and 
the decrement. The nominal values and errors used are given in Table ^ their description 
follows: 

D~£+ background fraction: The fractions of the two modes of D"^"*" background out of 
events in the signal peak (in the case of the semileptonic D~ decay channel it is done 
after subtraction of /other), y'^, and the fraction of Bd decays in each mode, 



and were calculated in Section |4.2.1| . Due to the dependance of and Y on the 



Bg production fraction fs, the uncertainty on the value of fs is a source of systematic 
uncertainty. 

B meson lifetimes, oscillation and production fractions: The world averages for 
the lifetime of the Bu, Bj and Bg mesons, for the B^ oscillation frequency and for the 
Bg production fraction were used. 

Oscillating combinatorial background fractions: The fractions of both types of oscil- 
lating combinatorial background out of the total combinatorial background for the 
various decay channels: Cq,Cq° '^,Cf and C°° were calculated in Section |6.3.1 . 



Other background fraction: The fraction of other types of background specific to semi- 
leptonic decays out of events in the signal peak: /other was estimated as described 



in Section 4.2.2 



Combinatorial background lifetime: The combinatorial background parameters 

/bg,iep' ^i,iep and Tbg^iep for semileptonic D" decays; and r^^had for hadronic D^" 
decays, were obtained from a fit to sideband events (Section |^. The statistical errors 
on the parameters from the sideband fit were used as systematic uncertainties. 

Mixing tag behavior in signal events: The uncertainties on the fraction of mixed events 
in each bin used to parameterize the deviation of the mixing tag from a true probability 
(as described in section |0 and shown in Figure Rl) include both statistical uncertainties 



from the limited Monte Carlo sample size, typically of order 0.025, and the following 
systematic effects: 
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Figure 9: Measured Bg oscillation amplitude as a function of Arris for this analysis. The 
error bars represent the la statistical uncertainties. The shaded bands show the ±1.645(T 
region, with and without including systematic effects. The values of Arrtg where the shaded 
region lies below the A = 1 line are excluded at the 95% confidence level. The dashed line is 
1.645a used to determine the experimental sensitivity, which is indicated by the circle. 
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• the effect of a one standard deviation variation in the Qsame offset, typically of 
order 0.004. 

• the effect of a one standard deviation variation in the Qopp offset, typically of 
order 0.005. 

• the effect of a one standard deviation variation in the fragmentation kaon tag's 
purity, typically of order 0.007. 

• the effect of a one standard deviation variation in the opposite lepton tag's purity, 
typically of order 0.006. 

• since the limited candidate sample size in data prevents us from showing that 
the simulation and data agree on the jet charge distributions with and without 
a fragmentation kaon tag, we take the entire effect of using those distributions 
instead of the single jet charge distribution as a systematic error, typically of the 
order of 0.03. 

The fitted uncertainties on the deviation parameters amix and /9mix were taken as the 
systematic uncertainties on the mixing tag. 

Mixing tag behavior in background: The fitted uncertainties on the deviations of the 
distributions of the mixing tag in combinatorial and backgrounds, relative 

to its behavior in signal (Section |6.3.2|) , were taken as additional systematic uncertain- 
ties. 

Decay length error correction: The fitted uncertainty on the decay length error correc- 
tion (Section |5.1|) was used as a systematic uncertainty. 

Decay time reconstruction bias: The fitted uncertainties on the residual bias in the de- 
cay time reconstruction (Section were used as systematic uncertainties. 

Detector resolution modelling: The resolution of the tracking detectors might affect the 
decay time reconstruction and the mixing tag. The simulated resolutions were degraded 
by 10 % relative to the values that optimally describe the data following the studies 
in Ijl^. The analysis was repeated and the mixing tag was found to be insensitive to 
this variation while the proper decay time resolution deteriorated by 5 %. This 5 % 
uncertainty on the proper decay time resolution was used as a bidirectional systematic 
uncertainty. 

The relative importance of the various systematic uncertainties, as a function of Arus, is 
shown in Table 0. For all ArUs values, the total systematic uncertainty is small compared 
to the statistical uncertainty. At low Arris the most important systematic contributions are 
from the uncertainties on the behavior of the mixing tag in the signal, while at high Arris 
the most important systematic contributions are from the uncertainties on the decay time 
reconstruction bias. 

10 Checks of the method 

10.1 Fitted decay length distribution and lifetime fit 

The form of the likelihood of a candidate decaying with a measured decay length /j has 
been stated in Section ^. By setting the various parameters to their nominal values and 
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input 
parameter 


nominal 
value 


Ci i \JL 


contribution to A.4. at Anis = 


ps~^ 


5 ps"^ 


10 ps"^ 


15 ps"^ 


Fractions of D background: 














0.143 


0.050 


1 n no 


+U.UO 


1 n no 

+u.uy 


1 n 1 1 
+U.ii 


n 


0.619 


0.045 


+U.UUii 


+u.uuiy 


1 n nn/i 
+U.UU4 


+U.UUO 




0.065 


0.035 


— U.UiU 


1 n non 


n nnn 

— u.uuy 


n noQ 
— U.Uoo 


fl 


0.470 


0.038 


1 n nni n 
+U.UUiy 


n nni n 
— U.UUiU 


n no7 
— U.Uz / 


1 n nni q 


Previously measured B meson properties: 










/. 


0.107 


0.014 


-0.030 


-0.028 


-0.03 


-0.030 




1.653 ps 


0.028 ps 


+0.0005 


-0.0004 


-0.013 


-0.0013 




1.548 ps 


0.032 ps 


+0.0021 


+0.0001 


-0.004 


+0.0008 




1.493 ps 


0.062 ps 


-0.004 


+0.004 


+0.016 


-0.04 




0.472 ps-^ 


0.017ps-^ 


+0.0009 


+0.0004 


-0.0002 


-0.0002 



Fractions of oscillating combinatorial background: 





0.163 


0.030 


-0.0033 


-0.004 


-0.015 


-0.018 


^no <l> 


0.270 


0.027 


-0.005 


+0.0032 


+0.0002 


+0.008 


cf 


0.163 


0.030 


-0.004 


-0.017 


-0.022 


-0.05 




0.049 


0.013 


-0.0029 


+0.007 


+0.010 


+0.04 



Fraction of other background: 
/other 0.135 



0.057 +0.006 +0.013 +0.010 +0.05 



Combinatorial background lifetime fit: 

J bg, lep 
bg, lep 



bg, lep 



bg, had 

Mixing tag behavior in signal events: 



0.92 


0.05 


+0.0020 


+0.0024 


-0.0016 


+0.011 


1.40 ps 


0.26 ps 


+0.0018 


-0.005 


-0.022 


-0.04 


0.7 ps 


0.5 ps 


-0.0013 


-0.007 


-0.012 


-0.024 


0.70 ps 


0.11 ps 


+0.014 


+0.011 


-0.04 


+0.06 





-0.006 


0.012 


-0.06 


-0.024 


+0.021 


Pmix. 


-0.051 


0.026 


-0.04 


-0.06 


-0.08 



Mixing tag distributions in background: 



combinatorial 


1 


0.18 


-0.017 


+0.0033 





Bu ^ B-e+ 


1 


0.14 


-0.007 


+0.0023 


-0.0023 



+0.005 
-0.06 

-0.04 
-0.011 



Decay length resolution: 



fitted correction 
detector modelling 


1.478 
1 


0.033 
0.05 


+0.0009 
+0.0012 


+0.016 
+0.035 


-0.005 
-0.010 


+0.06 
+0.14 


Decay time bias: 

''bias 

.slope 

''bias 


-0.033 ps 
1.29 


0.005 ps 
0.16 


+0.0008 
+0.0010 


-0.006 
-0.005 


+0.031 
+0.009 


+0.33 
+0.07 


systematic uncertainty 
statistical uncertainty 


0.12 
0.28 


0.11 
0.64 


0.15 
1.2 


0.41 
3.0 



Table 2: Input parameters for the fit and their contributions to the systematic errors, 
as described in the teoct. Systematic errors of less than 5 • 10~^ are reported as zero. For 
comparison the total systematic errors and the statistical errors are given as well, in the last 
rows. 
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the amplitude to zero, we acquired a prediction for the distribution of k. This prediction 



compared well with the actual measured values of /j, as shown in Figure \L0 . 

The likelihood can also be used to measure the Bg lifetime by ignoring the mixing tag and 
fitting the Bg lifetime to the data. The resulting lifetime is 1.57=b0.17ps, which is consistent 
with the world average value of 1.54 ± 0.07 ps. 

10.2 Oscillation fit on simulated events 

A likelihood fit for Arris was performed on a simulated Monte Carlo event sample having 
the same statistics and estimated composition as the data, with an oscillation at a true 
frequency of either Amg = 2.0 ps^^ or Amg = 3.0 ps^^. The fit yielded Arris = l-98lo;J| ps~^ 
and Arris = 3.16lg;| ps respectively, in agreement with the input values. Performing an 
amplitude fit (ignoring systematic uncertainties) on the same Monte Carlo events yielded 



the results shown in Figure |TT[ As expected, the amplitude is consistent with 1 at the true 
value of Arris. 

The sensitivity of the analysis is defined as the expected highest oscillation frequency 
excluded at the 95% confidence level, given that the true Arris is infinitely high. Given an 
infinitely high Arris the expectation value of the amplitude at all fitted values of Arris is 
zero 0. This allows us to evaluate the sensitivity as the frequency at which the result- 
ing 1.645(7 line rises above an amplitude of one, as shown in Figure ^. The experimental 
sensitivity of this analysis is 4.1 ps~^. 



11 Conclusion 

A sample of Bg decays obtained using Dg"£+ combinations was used to study Bg oscilla- 
tion. The estimated sensitivity of the analysis is 4.1 ps~^. The resulting, stand-alone, lower 
limit from our analysis is significantly lower, at 1.0 ps~^. This limit is not competitive with 
existing limits. However, this does not diminish the contribution of this analysis to the 
world combined measurement, which is at the rapid oscillation region. The previous OPAL 
measurement was an inclusive measurement, and so has only a negligable statistical corre- 
lation with our measurement. Its sensitivity was 6.7ps~^ and it set a lower limit of 5.1 ps~^. 



Combining it with our results we get the combined measurement shown in Figure |1^, a 
sensitivity of 8.0 ps~^ and a limit of 5.1 ps^^. 
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Figure 10: Comparison of expectation from fit results and actual decay length distributions. 
The fit results are shown by the curves, with individual components shaded, while the decay 
lengths measured in the data are shown by the error bars. The sideband region was used in 
the combinatorial background lifetime fit. 
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Figure 11: Measured Bg oscillation amplitude as a function of Amg in simulated events 
with a true Amg o/2.0ps~^ in the upper plot, and 3.0 ps~^ in the lower plot. The error 
bars represent the la statistical uncertainties. The shaded bands show the ±1.645(j region 
(systematic effects not included). The arrows point at the simulated value of Arris. 
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